Composite reweighting with Imaginary Chemical Potentials in SU(3) 
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We review the overlap pathology of the Glasgow reweighting method for finite density QCD, and discuss the 
sampling bias that effects the determination of the ensemble-averaged fugacity polynomial expansion coefficients 
that form the Grand Canonical Partition function. The expectation of the difference in free energies between 
canonical partition functions generated with different measures is presented as an indicator of a systematic quark 
number dependent biasing in the reweighting approach. The advantages of building up an unbiased polyno- 
mial expansion for the Grand Canonical Partition function through a series of parallel ensembles generated by 
reweighting with imaginary chemical potentials are then contrasted with those of addressing the overlap pathology 
through a secondary reweighting. 



1. Introduction 

Interest in the numerical determination of 
BCS-type instabilities at the Fermi surface in 
QCD at high density has led to a recent resur- 
gence of interest in the analytic continuation of 
Monte Carlo analyses and reweighting . For 
QCD the SU(3) lattice action |^ becomes com- 
plex at nonzero chemical potential /i 7^ 0, and so 
the probabilistic importance sampling steps of a 
Markov chain are undefined prohibiting standard 
Monte Carlo techniques. The Glasgow method 
has given a useful means of quantifying the suc- 
cess of these approaches, where the /i dependence 
of the Grand Canonical Partition function r2(/i) 
is made semi-analytic through polynomial expan- 
sion in fugacity z = e'^/^ 

Although the complex action issue is negated 
in generating a lattice ensemble in an accessible 
regime of the parameter space (eg. fi = for 
SU(3)), the analysis is pathological. An unphys- 
ical onset transition in the quark number den- 
sity Jo 3X fx = ^TOtt is thought be a feature of 
quenched (or isolated configurations of) QCD at 
finite density, but in addition exponentially high 
statistics are required to generate physical results 
from dynamically generated ensembles. It is this 
second overlap pathology which also persists for 
QCD-like models with the Glasgow scheme in 
which the Goldstone modes of the theories are ei- 
ther degenerate or considerably heavier than \m-,^ 
I^l 1^ , which hints further at a statistical origin. 

By evaluating the ratio of ensemble-averaged 



terms in these QCD-like models (for which the 
fermion determinant is positive definite) we were 
able to study the comparative effect of reweight- 
ing with the Glasgow method for ensembles gen- 
erated at different values of chemical potential /io, 
and in so doing develop a new reweighting ap- 
proach for QCD. 

2. Glasgov^f Method 

The characteristic polynomial in fugacity is 
formed through the definition of the propagator 
matrix P |^ from gauge links in the spatial G and 
time directions V . 
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Since V is an overall factor of P the order of the 
expansion can be reduced by exploited the uni- 
tary symmetry Z„j defined by e^'^*^^"* (where j 
is integer). This relates the expansion terms to 
the GCPF through the Canonical Partition 
functions for n quarks, represented by the suit- 
ably normalised ensemble-averaged expansion co- 
efficients in a reweighted measurement. 
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In addition to possible autocorrelations 
through the numerical evaluation of the ensemble 
there is also an inherent covariance between the 
numerator and denominator of a reweighted mea- 
surement which becomes exponentially dominant 
for significantly different ratios of weightings ||] . 
The Cauchy distribution of the ratio population 
formed by two independent gaussian distribu- 
tions gives a simplistic picture of this behaviour. 
The joint density function integral for this ratio 
population diverges and the higher moments are 
thus undefined. Only for identical numerator and 
denominator populations does the joint density 
function then reduce to a gaussian, following 
suitable substitution. 

For SU(3) given the accessible region of the 
T — fi phase plane for importance sampling meth- 
ods, the ratio of weightings given in Eqn.8 gen- 
erally does differs from one when the expansion 
is reweighted to the vicinity of the critical line. 
Kurtosis in the measurement from the weighted 
average relative formed with the ratio of weight- 
ings is of particular concern for this scheme since 
the expansion coefficients vary over several hun- 
dred orders of magnitude. In SU(3) for a 4^ 
volume typically at /3 = 5.04, c„=o ^ £^ and 
and the kurtosis is thus uneven 
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across the expansion. For the limited range of 
the suitably normalised expansion terms of order 
one (from an ensemble generated at some /ii), the 
numerator and denominator populations become 
equivalent and the averaging is then independent 
of the reweighting evaluation at /Xq. Since these 
expansion coefficients are in general complex inte- 
grating over all gauge configurations cancellation 
occurs of all but the triality zero terms (rnod{j, n) 
j — 1, 2), but with a significant real- valued kurto- 
sis this cancellation becomes inexact. The phase 



of the ratio of complex weightings is of less conse- 
quence given this Z{3) tunneling and the equality 
of the transition probability for conjugate gauge 
configurations. 

The related effect was noticed in the evalu- 
ation of Q{ih') by reweighting using dynamical 
ensembles generated at various values of ivg for 
the Hubbard model To address the kurto- 

sis with complex weightings for QCD we could 
envisage recovering each fl n through the Fourier 
transform of using a reweighted fugacity 

expansion with an ensemble generated at some 
ivo- With this scheme for 17 (iz/) to be reliably 
sampled through reweighting we would have to 
evaluate ensembles for a series of values of ivo, 
but the only non-vanishing term of the fugacity 
expansion from the Fourier transform is cjv, and 
the identification of a unique value of i^o for which 
the sampling of c„ is effective is therefore explicit. 
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However, rather than evaluate a costly series of 
ensembles generated for a range of values of ivo to 
ensure the analyticity of the fugacity polynomial 
we can make explicit measurements for the nu- 
merator and denominator to address the kurtosis 
of complex weightings. The eigenvalue problem 
to identify /iti for a given c„ is first solved, and 
then the numerator evaluated with these normali- 
sations for each term of the expansion. The prod- 
uct of the weighting ratios for each consecutive 
value of /ii is then evaluated, thus ensuring the 
denominator evaluation is also of order one. 

Our results for 10,000 configurations of a 4^ 
lattice with four staggered quark flavours at /3 = 
5.04 and m = 0.1 are presented in Fig.l, with 
the Lee- Yang zeros of the GCPF evaluated in the 
complex-yit plane. Our improvement scheme re- 
veals a rich finite density phase structure, with 
two clear transition points (where the zeros ap- 
proach the real axis). The first is in agreement 
with In) ( Rc(/3o)=4.98(9) ), and the second also 
identified in p2[ appears from our preliminary 
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Figure 1. Lee- Yang zeros for QCD evaluated in 
th e complex-/^ plane for a 4^ lattice aX (3 — 5.04 . 



scaling analysis to be associated with lattice sat- 
uration. Notably, the critical feature at /i ~ 0.4 
appears to correspond to imTr for this bare mass 



13 1, and it will be interesting to repeat the anal- 
ysis at smaller values of m. 

3. Summary 

Each analytic continuation scheme reliant on 
reweighting can be pictured as a Taylor expansion 
about the zero density ensemble, and through the 
Glasgow scheme we can identify the range of va- 
lidity of ensemble-averaging for complex weight- 
ings. We have addressed the inherent kurto- 
sis of measurement by redefining our ensemble- 
averaging and can now more reliably evaluate the 
Grand Canonical Partition function, ensuring an- 
alyticity in the fugacity variable. Using in addi- 
tion reweighting simultaneously in (3 we could also 
now reliably reweight to the entire T— /j, plane (up 
to lattice saturation for finite volumes) given the 
independence on Monte Carlo measure we have 
established. 

The Lee- Yang zeros we have evaluated from our 
new scheme appear to approach the real axis in 
manner accordant with the finite volume scaling 
circle theorem of 10. However, it remains to be 
seen if we can make useful comments on diquark 
condensation at finite density, and evaluate the 
critical exponents of this transition following ap- 
proaches suggested in p5||lql through the finite 



volume scaling of the edge singularity revealed 
with this analysis. 
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